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Abstract. In this article, we prove the existence of Arnold diffusion for an in- 
teresting specific system — discrete nonlinear Schrodinger equation. The proof 
is for the 5-dimensional case with or without resonance. In higher dimensions, 
the problem is open. Progresses are made by establishing a complete set of 
Melnikov-Arnold integrals in higher and infinite dimensions. The openness 
lies at the concrete computation of these Melnikov-Arnold integrals. New ma- 
chineries introduced here into the topic of Arnold diffusion are the Darboux 
transformation and isospectral theory of integrable systems. 
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1. Introduction 

For a simple example posed by V. I. Arnold [5], the existence of the so-called 
Arnold diffusion has been proved (see e.g. |3] @] |14|h The argument involves two 
parts: A calculation of the Melnikov-Arnold integrals |3] and a transversal inter- 
section argument j4] supported by a A- lemma |14) . Other arguments of variational 
type were also developed ED El 0- 

Nevertheless, the theory of Arnold diffusion is far from complete |13] [5] |28) 
EE! EH] © M E2] EH US! EB The main challenge is dealing with high di- 
mensional specific systems of interest in applications. When the dimensions of the 
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KAM tori are large, more Melnikov- Arnold integrals are needed to establish Arnold 
diffusion. Calculating these integrals is a dauting task if not impossible, even with 
computers. In an infinite dimensional phase space, even the dimensions of the KAM 
tori become a challenging issue. Are there infinite dimensional KAM tori ? in what 
form ? In the classical setting of a Banach space with angles-momenta coordinates, 
the challenge is how to deal with the perturbation v.s. the decay of the sequence 
of momenta. It is a very interesting problem. 

The aim of the current article is to draw attention to two canonical systems of 
mathematical physics: The discrete nonlinear Schrodinger equation (DNLS) and its 
continuous version - the nonlinear Schrodinger equation (NLS). DNLS and NLS are 
integrable systems that describe many different phenomena in physics [2]. DNLS 
is an integrable finite difference discretization of NLS. An interesting fact about 
DNLS is that one can choose and change the dimensions of the phase space by 
selecting the number of particles in the discretization. For a two particle case, 
under periodic Hamiltonian perturbations, the resulting system is 5-dimensional, 
for which the existence of Arnold diffusion will be proved here. For a three particle 
case, the system is 7-dimensional, and will be a good testing ground for a higher 
dimensional theory. 

The integrable theory offers the missing link from low dimension to high di- 
mension via two powerful and beautiful machineries: Darboux transformation and 
isospectral theory. Darboux transformation generates explicit expressions of sepa- 
ratrices, while isospectral theory produces all the Melnikov vectors. Together they 
provide a complete set of Melnikov- Arnold integrals with elegant universal formu- 
lae. This is the case for both DNLS and NLS dU HU US] EO] |2l]. On the 
other hand, specific calculation of these integrals is the challenge. For the purpose 
of proving the existence of chaos, often one Melnikov integral is enough and eas- 
ily computable [19] |21| |22| [25] [26) . The reason is that one can utilize locally 
invariant center manifolds instead of KAM tori. For the purpose of Arnold diffu- 
sion, local invariance (i.e. orbits can only enter or leave the submanifold through 
its boundaries) is not enough. This is due to the second part of the argument 
for Arnold diffusion mentioned above. To establish a A-lemma, one needs the tori 
to be invariant (not locally invariant). For a different application of A-lemma in 

(1) . establishing shadowing lemma in infinite dimensional autonomous systems, 

(2) . proving the existence of homoclinic tubes and heteroclinically tubular cycles, 

(3) . proving the existence of tubular chaos, (4). proving the existence of chaos 
cascade; we refer the readers to |23| |24| |25j |26| . 

The article is organized as follows: Section 2 deals with isospectral theory and 
Darboux transformation for both DNLS and NLS. Section 3 deals with Arnold 
diffusion for a 5-dimensional perturbed DNLS. 

2. Isospectral Theory and Darboux Transformation 

In this section, we are going to present the isospectral theory and Darboux 
transformation for both the discrete nonlinear Schrodinger equation (DNLS) and 
the nonlinear Schrodinger equation (NLS). 

2.1. Discrete Nonlinear Schrodinger Equation. Consider the discrete 
nonlinear Schrodinger equation (DNLS), 

(2.1) iq n = To [g n +l - 2 <7« + Qn-l] + l<7n| 2 (<?n+i + Qn-i) - 2u 2 q n , 
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where q n 's are complex- valued, i = y— T is the imaginary unit, w is a positive 
parameter, and q n satisfies the periodic boundary condition and even constraint, 



in , 



(2.2) q n+N = q n , q- n = <?„ 

where TV is a positive integer N > 3 and h = l/N. The DNLS is a 2(M + 1)- 
dimensional system, where M = N/2 (N even) and M = (N - l)/2 (JV odd). The 
DNLS can be rewritten in the Hamiltonian form 

(2.3) iq n = p n dH /dq^ , 
where p„ = 1 + h 2 \q n \ 2 and 

1 N ~ 1 (_ 2 1 

^0 = 7^^ { 1n{ln+l + ?n-l) - ^(l + W 2 /l 2 )lnp„ > . 

n=0 

The phase space is defined as 



S = \ q = (q,q) 



= (g , 9i, -,qN-i), QN-n =Qn {1 < n < N - 1) 

In S (viewed as a vector space over the real numbers), we define the inner product, 
for any two points q + and q ~, as follows: 

N—i 

(q + , q ") = {qtqn + qnQn) ■ 

And the norm of q is defined as \\q || 2 = (q,q). 

Remark 2.1. In the expression of i? , both X^^To 1 (<7n+i + 3n-i)] an d ^ = 
W ~52n=o m P« are constants of motion too. / will be used later to establish Arnold 
diffusion in the non-resonant case. Also the constant of motion D given by D 2 — 
Y[n=o P n wm pl a y an important role in the isospectral theory. In the continuum 
limit (i.e. h — ► 0), the Hamiltonian Hq has a limit in the manner: hHo —* H c , 
where H c is the Hamiltonian for NLS, H c = — L [\q x \ 2 + 2uj 2 \q\ 2 — |g| 4 ]c?x. Also as 
h->0, D->1. 

2.2. Isospectral Theory of DNLS. For more details on the topic of this 
subsection, see |18j . DNLS has the Lax pair [J 

(2.4) ip n+1 = L n <p n , 

(2.5) (f> n = B n (fi n , 



where 



z ihq„ 
ihqn~ l/z 

i { 1 — z 2 + 2i\h — h 2 q n q n -i + uj 2 h 2 —izhq n + (1 / z)ihq n -i 



h 2 V -i^hq^Zi + (l/z)ihq~n~ p - 1 + 2iXh + h 2 q^q n ^i - u> 2 h 2 

and where z — exp(iA/i). Compatibility of the over determined system 12.412.5(1 
gives the "Lax representation" 

L n = B n +iL n — L n B n 

of the DNLS (|2.1|) . Focusing our attention upon the discrete spatial part (|2.4I) of 
the Lax pair, let M n = M n (z,q) be the 2x2 fundamental matrix solution such 
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that Mo is the 2x2 identity matrix. The Floquet discriminant A : C X S — ► C is 
defined by 

A = — trace MN(z,q) , 

where £> 2 = rin^Q 1 / 9 "- The isospectral theory starts from the fact that for any 
z G C, A(z, g) is a constant of motion of DNLS. An easy way to understand this 
is that as q n evolves in time according to DNLS, the parameter z in the Lax pair 
does not change. A(z, q) is a meromorphic function in z of degree (+N, —N), and 
provides (M + 1) functionally independent constants of motion, where M = iV/2 
(TV even), M = (N — l)/2 (JV odd). There are many ways to generate (M + 1) 
functionally independent constants of motion from A. The approach which proved 
to be most convenient is by employing all the critical points z c of A(z,q): 

<z c ,q) 



dz 



. 



Definition 2.2. We define a sequence of (M 
follows 



1) constants of motion Fj as 



(2.6) 



Fj (q) = A(zUq), q), (J = 1, • • • , M + 1) 



There is a good description on the locations of these critical points z| in the 
NLS setting |19j . These Fj's can be used to build a complete set of Melnikov- 
Arnold integrals for the Arnold diffusion purpose. The Melnikov vectors are given 
by the gradients of these Fj 's. 



Theorem 2.3. Let z c Aq) be a simple critical point, then 



(2.7) 



dFj 
dq n 

dFj 



I ^^ 2) +4 + ' 2u( -' 2i 



2W n 



+1 



+1 



V 



^(+,1) -(-,1) 



wher e tp± = (V'i ± ' 1 \ -0i ±,2) ) T ; (T: 
pair \2.412.5)) at (Zj,q) such that 



transpose), are two Bloch solutions of the Lax 



where ^ are periodic in n with period N , ( is a complex constant, and W n = 
det (V'ni'^n) * s the Wronskian. 

For a perturbation Ho + eHi of the Hamiltonian Ho, a complete set of Melnikov- 
Arnold integrals is then given by 



,00 N-l 
n =Q 



BR 



dq n defc 



Pn 



dt 



, M + 1) 



which arc evaluated on the unstable manifolds of tori that persist into KAM tori, 
where is given by l|2.7|l . Expression 12.7|l is a universal expression for all the 
Melnikov vectors. The challenge is of course how to compute these Ij's. The 
difficulty lies at how to obtain expressions of the orbit q n and the corresponding 
"0^. By utilizing Darboux transformations, these expressions can be obtained even 
by hand in some cases. 
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2.3. Darboux Transformation of DNLS. For more details on the topic of 
this subsection, see |18| . Expressions of unstable manifolds of tori can be gener- 
ated via Darboux transformations. For DNLS, such a Darboux transformation was 
established in |17| . 

Let q n be any solution of DNLS. Pick z = z d (\z d \ ^ 1) at which the linear op- 
erator L n in the Lax pair H2.4|) - H2.5(l has two linearly independent periodic solutions 
(or anti-periodic solutions) <p^ in the sense: 



(or <t>t +N = -Dtfc) , 



where D is defined above. Let 



i>nl,4>n2) T = C 4 



where c + and c are complex parameters. Define the matrix r„ by 

r„ = 



z+(l/z)a n b n 

C n -l/z + zdn 



where 



z d2 A n 



\0n2\ 2 + \z d \ 2 ttnl\ 2 



d„ 
b,, 



1 



z d A 
z d \ A - 1 



6„ 2 r + kTi<M 



z d A r 

\z d \ 4 ~ I — 



Unl9n2, 



\z d \ 2 A., 
1 



Unl<Pn2, 



nlV + \z d \ Z \<j)n2\ 



From these formulae, we see that 

tin — d nj b n — C n - 

Theorem 2.4 (Darboux Transformation [T7]). Define Q„ and \P n by 

Qn "T^n+l ®n+lQn 3 TriV^n 3 

h 

where ip n solves the Lax pair \2.4\) - \2.5)) at (q ni z). Then Q n is also a solution of 
DNLS, and ^ n solves the Lax pair \H.4^ - ^-5\) at (Q n ,z). 

In principle, by choosing q n to be orbits on the tori, one can generate Q n to be 
orbits on the unstable manifolds of the tori. In special cases, Q n can be calculated 
by hands. 

Next we present an example. Define the 2-dimensional invariant plane 

(2.8) IT = {if G <S I q n = q, Vn} . 

On II, the solutions of DNLS are given by the periodic orbits (1-tori) 

(2.9) q n = q c ,Vn; q c = aexp {~i[2(a 2 - uj 2 )t - 7]} , 
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where a and 7 are real constants. We choose the amplitude a in the following range 
so that the unstable direction of q c is 1-dimensional 



TV tan — < a < N tan — 
N N 



when N > 3 



(2.10) 



3 tan — < a < 00 



when N = 3 . 



Increasing the unstable dimensions of g c amounts to iterations of the Darboux trans- 
formation which are still doable by hands, and does not add difficulty substantially 
in the Arnold diffusion problem. To apply the Darboux transformation we choose 



z = Vpcos— + ^/pcos 2 — - 1 , 

where p = 1 + h 2 a 2 . This z d is also a critical point of A. We label it by z\ = z d . 
Direct calculation leads to the following formulae 



(2.11) 

(2.12) 
where 



Q71 



r/A„ 



dF\ 
dq n 

dF\ 



KiKn}- 1 scch[2 M t + 2p] I £ 2 



r = 

A„ = 

K = 

K n = 



1 — cos 2ip — i sin 2ip tanh[2/it + 2p] , 

1 ± cost/? [cos /3] _1 sech[2^i + 2p] cos[2n/3] , 

-2N(1 - z 4 )[8ap 3/2 z 2 ]-Vpcos 2 /3- 1 , 

cos (3 ± cos ip sech[2fit + 2p] cos[2(n — l)/3] 



cos (3 ± cos ip sech[2/^i + 2p] cos[2(n + l)/3] 
cos (3 sech [2/zt + 2p] ± (cos </? 



A„ 



— « sintp tanh[2/xi + 2p\) cos[2n/3] 

cos /? sech [2/rf + 2p] ± (cos <p 
+i sin tanh[2^it + 2p\) cos[2n/3] 



i20(t) 



i20(i) 



P = n/N, [i = 2/^ 2 Vpsin/Vpcos 2 /? - 1 , 

p = 1 + /i 2 a 2 , h=l/N, z = yfp cos (3 + \fp cos 2 (3 



1 , 



0(t) = (a 2 - w 2 )t - 7/2 , ftae iv = Vpcos 2 /?- l + i^/psin/3 . 

and p is a real parameter. One can easily see that Q n represents homoclinic orbits 
asymptotic to the periodic orbits q c : As t — > ±00, 

i(7T±2¥>) 
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The union 

IJ Qn 

7e[0,27r] 

represents the 2-dimensional unstable (=stable) manifold of the 1-torus (|2.9|l . When 
N > 3, the 1-torus also has a center manifold of codimension 2. 

2.4. Nonlinear Schrodinger Equation. Consider the nonlinear Schrodinger 
equation (NLS), 

(2.13) iq t = q xx + 2[\q\ 2 -uj 2 ]q , 

where q = q(t, x) is a complex- valued function of the two real variables t and x, 
t represents time, and x represents space. q(i, x) is subject to periodic boundary 
condition of period 2-7T, and even constraint, i.e. 

q(t, x + 2n) = q(t, x) , q(t, —x) — q(t, x) . 

a; is a positive parameter. The DNLS i|2.1[l is an integrable finite difference dis- 
cretization of the NLS. The NLS can be rewritten in the Hamiltonian form 

(2.14) iq = dH /dq , 
where 

H = - r [\q x \ 2 + 2oj 2 \q\ 2 -\q\ 4 ]dx . 
Jo 



10 

The phase space is defined as 

S = lq= (q,q) 

where H k is the Sobolev space of periodic and even functions. 



qeH k , (fe>l)| 



2.5. Isospectral Theory of NLS. For more details on the topic of this sub- 
section, see |19| . NLS has the Lax pair 

(2.15) ^ = Ui/> , 

(2.16) Vt = V^j, 



where 



A q 

q -X 



V 



U = i 

2A 2 - \q\ 2 + lu 2 2\q-iq x 



2\q + iq x - -2A 2 + \q\ 2 - uj 2 

Focusing our attention on the spatial part 12.15fl of the Lax pair (|2. 1512. 1611 . we can 
define the fundamental matrix solution M(x) such that M(0) is the 2x2 identity 
matrix. Then the Floquet discriminant A is defined as 

A = trace M (2tt) . 

The isospectral theory starts from the fact that for any A S C, A(A, q) is a constant 
of motion of NLS. A = A(A, q) is an entire function in both A and q. A provides 
enough functionally independent constants of motion to make the NLS H2.13|l in- 
tegrable in the classical Liouville sense. There are many ways to generate these 
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functionally independent constants of motion from A. The approach which proved 
to be most convenient is by employing all the critical points A c of A(A, q): 

For each q, there is a sequence of critical points {A^} whose locations can be esti- 
mated nnj. 

Definition 2.5. We define a sequence of constants of motion Fj as follows 

F j (q)=A(X%q),q) . 

These F/s can be used to build a complete set of Melnikov- Arnold integrals 
for the Arnold diffusion purpose. The Melnikov vectors are given by the gradients 
of these Fj's. 

Theorem 2.6. Let A^(g) be a simple critical point, then 



(2.17) 



dq _ . VA2-4 / ifj+lj)- 

9Fj ~ l w(ip+,i>-) { -v>i>r 



poo p2ir 


" dFj dHi 




/ / 2 Im 




dxdt 


l-oo JO 


dq dq _ 



where ?/' ± = (ip^, 4'2 : ) T , (T : transpose), are two Bloch solutions of the Lax pair 
JpJh}) - Jfilb}) at (X c p q) such that 

where a is a complex constant and tp^ (x) are periodic in x with period 2n, 

is the Wronskian, and A is evaluated at A = A|. 

For a perturbation Ho + eHi of the Hamiltonian Ho, a complete set of Melnikov- 
Arnold integrals is then given by 

h = 

which are evaluated on the unstable manifolds of tori that persist into KAM tori, 
where -g*- is given by (|2.17(l . Expression Ij2.17|l is a universal expression for all 
the Melnikov vectors. The challenge is of course how to compute these I/s. The 
difficulty lies at how to obtain expressions of the orbit q and the corresponding tp^ . 
By utilizing Darboux transformations, these expressions can be obtained even by 
hand in some cases. Also as mentioned in the introduction, there is the issue of 
elusiveness of infinite dimensional KAM tori. As j — > oo, The magnitude of Ij v.s. 
the size of the perturbation e is another tricky problem. 

2.6. Darboux Transformation of NLS. For more details on the topic of 
this subsection, see 119) . Expressions of unstable manifolds of tori can be generated 
via Darboux transformations. For NLS, such a Darboux transformation was also 
known. 

Let q be any solution of NLS. Pick A = v (complex constant) at which (|2.15|) 
has two linearly independent periodic solutions (or anti-periodic solutions) ± in 
the sense: 

± (x + 2tt) = ^{x) , (or ± (a; + 27r) = -^(a:)) . 
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Let 

= c + <j) + + c_0~ , 
where c+ and c_ are complex parameters. Define the matrix G by 

where 

p _ A 01 -02 \ 
V 02 01 / 

Theorem 2.7 (Darboux Transformation). Define Q and "J 6y 

Q = g + 2(y - p) |^ 1 -H^ ' * = G ^' 

where ip solves the Lax pair y2.15)) - \2.1b)) at (g, A). TTien Q * s a ^ so a solution of 
NLS, and ^ solves the Lax pair \2.15)) - \2.1b)) at (Q, A). 

In principle, by choosing q to be orbits on the tori, one can generate Q to be 
orbits on the unstable manifolds of the tori. In special cases, Q can be calculated 
by hands. 

Next we present an example. Define the 2-dimensional invariant plane 

n = {qes\ d x q = o} . 

On II, the solutions of NLS are given by the same periodic orbits (1-tori) as in (|2.9|) 



(2.18) q c = aex-p{-i{2(a 2 -uj 2 )t-j}} , 

where a and 7 are real constants. We choose the amplitude a in the range a € 
(1/2, 1) so that the unstable direction of q c is 1-dimensional. Increasing the unstable 
dimensions of q c amounts to iterations of the Darboux transformation which are 
still doable by hands |21| . and does not add difficulty substantially in the Arnold 
diffusion problem. To apply the Darboux transformation we choose 



a 



Va? - 1/4 



This v is also a critical point of A. We label it by z\ = v. Direct calculation leads 
to 

r 1 _1 

Q = q c 1 ± sin -do sechr cos x 



cos 2i? — i sin 2i? tanh r =F sin i9 sechr cos x 



c)Fi 



dq 



where 



iti 

"2 



4 (\uiY + \u 2 YY 



cosh — cos z — i sinh — sin z , 
2 2 

7- T 

— sinh — cos(z — $ ) + i cosh — sin(z — ?9 ) , 



q c ui 2 
-q c U2 



2 



- + z/, T = 2crt-p, z = ac/2 + tf /2 =f tt/4 
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and p is a real parameter. As t — > ±00, 

Q - q c e^° . 

The union 

U Q 

7e[0,27r] 

represents the 2-dimensional unstable (=stable) manifold of the 1-torus l|2.18[l . The 
1-torus also has a center manifold of codimcnsion 2. 

3. Arnold Diffusion 

To establish the existence of Arnold diffusion, one needs three ingredients: 
(1). Melnikov-Arnold integrals, (2). A A-lemma, (3). A transversal intersection 
argument. Melnikov-Arnold integrals have been studied above. Next we discuss 
the other two ingredients. 

Lemma 3.1 (The A-lemma of Fontich-Martin [l4j). Let F be a C 2 diffeomor- 
phism in R", T be a C 2 invariant torus in M™. The dynamics on T is quasi-periodic. 
T has C 2 unstable, stable and center manifolds (W u , W s , W c ). Let V be a C 1 man- 
ifold intersecting transver sally W s at a point. Then 

W u C |J F m (T) . 

m>0 

Remark 3.2. Like every other A-lemma, the claim is very intuitive, but the 
proof is always delicate. The proof in |14j takes about eight pages. For a quick 
glance of the basic idea, see e.g. |23j . It turns out to be crucial to use Fenichel's 
fiber coordinates. With respect to base points, Fenichel fibers drop one degree 
of smoothness. That is why C 2 smoothness is required in the lemma to obtain 
C 1 families of C 2 unstable and stable Fenichel fibers. Using the Fenichel's fiber 
coordinates, W u and W s are rectified, i.e. they coincide with their tangent bundles. 
This makes the estimate a lot easier. Fenichel fibers also make it easier to track 
orbits inside W u and W s via the fiber base points in T. The main argument is 
to track the tangent space of a submanifold S of T starting from the intersection 
point of r and W s . After enough iterations of F, one can obtain some estimate of 
closeness to W u . The novelty of |14| is that they also track the tangent space at 
every point in S and off W s . This is necessary in order to obtain the claim of the 
lemma. The claim is proved by showing that any neighborhood of any point on 
W w has a nonempty intersection with F m (S) for some m. The claim of the lemma 
should also be true in proper infinite dimensional settings. 

Definition 3.3 (Transition Chain). A finite or infinite sequence of tori {Tj} 
forms a transition chain if W s (Tj) intersects transversally W u (Tj+i) at some point, 
for all j, and dynamics on Tj is quasi-periodic. 

Lemma 3.4 (Arnold Let {Tj} (I < j < N) be a finite transition chain. 
Then an arbitrary neighborhood of an arbitrary point in W u (Ti) is connected to an 
arbitrary neighborhood of an arbitrary point in W s (Tn) by an orbit. 



ARNOLD DIFFUSION 



245 



Proof. Let fli be an arbitrary neighborhood of an arbitrary point u% <E 
W u (Ti). Let B Ul (n) C fli be a closed ball of radius n > centered at u\. 
Then using the A-lemma l3.ll one can find 

(3.1) B UN {r N ) C S UJV _ 1 (r A r_i) ■ • • C B U2 (r 2 ) C B U1 {t\ ) 

such that 

Uj eW u (Tj), rj >o, (l<j<N). 
Indeed, by the A-lemma l3.ll 

B Ul {rx)r\W u {T 2 )^% 
where B Ul (n ) is the open ball. Then one can find 

u 2 G B„,(ri) n W U (T 2 ) , and 5„ 2 (r 2 ) C B Ul (n) . 
Again by the A-lemma l3.ll 

B U2 (r 2 )nW u (T 3 )^(b , 

and one can find 

u 3 G B„ 2 (r 2 ) n W U (T 3 ) , and S tt3 (r 3 ) C B„ 2 (r 2 ) , 

and so on. Let £In be an arbitrary neighborhood of an arbitrary point in W s (Tn). 
Inside fijy, one can find a C 1 submanifold intersecting transversally W s (Tn) at the 
point. Again by the A-lemma l3.ll 

B UN (r N ) n F m {Q N ) ^ for some to . 

Thus r2i and fljy are connected by an orbit. □ 

Remark 3.5. It is easy to see that around the connecting orbit, there is in fact 
a connecting flow tube [57 ] [5i | [25 ] [25 ] 

|J F m (£>) for some M 

0<m<M 

where D is a neighborhood. When TV = oo, relation l|3.1[) leads to a point in the 
intersection 

oo 

Starting from u, one obtains a connecting orbit. In a Banach space setting, if 
N = oo, then one can choose rj+i < f° r an y 3 m 3 - II) - Choosing an arbitrary 
point Vj in B Uj (rj), one gets a Cauchy sequence Thus 

oo 

lim vj = v€ H B u (r 3 -) . 
Starting from u, one still obtains a connecting orbit. 
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3.1. Arnold Diffusion of DNLS (N = 3, Non-resonant Case). In this 
subsection, we prove the existence of Arnold diffusion for a perturbed DNLS when 
N — 3, which is a 5-dimensional system. For arbitrary N, one can find large enough 
annular region inside the invariant plane II l|2.8|) . which is normally hyperbolic, 
for which the current proof can be easily applied. The point is that increasing 
unstable and stable dimensions does not pose substantial computational difficulty 
to establishing Arnold diffusion, while increasing the dimensions of tori does. We 
will study here the case that there is no resonance (uj = 0) inside the invariant 
plane II Q2.8p . The resonant case (u> ^ 0) will be studied in next subsection. 

Consider the following perturbation of the DNLS l|2.I() 

H = H a + efli , 



where 



Hi = a sin t 



N-l 

E 

n=0 



q n -i 



N-l 

E 



Qn-l 



q n -i 



where a is a real parameter. Under this perturbation, dynamics inside LI is un- 
changed, fl consists of periodic orbits forming concentric circles (|2.9ll [Figure [I]. 
We are interested in the following normally hyperbolic annular region inside II 



A 



q, Vn, 



3tan| < \q\ < Bj 



where B is an arbitrary large constant. Denote by {T u ' s {q) ■ q G A} the C 1 
families of C 2 one dimensional unstable and stable Fenichel fibers with base points 
in A TU] such that for any <f» S F u {q) or £ F s {q), ( q& A), 

WF^q*) - 511 < Ce Kt \\q» - q\\ , V< £ (-oo,0] , 

or 

||f(&) - q\\ < Ce- Kt \\ql - q\\ , V* G [0, oo) , 
where F* is the evolution operator of the perturbed DNLS, k and C are some 
positive constants. The Fenichel fibers are C 1 in e G [0, eo) for some eo > 0. It 
turns out that the constant of motion of DNLS (|2.1|l 

1 N ^ 
I= E ln ^« 

Tl = 

and Fi [cf: Ij2.6|l and l|2.12|l ] are the best choices to build the two Melnikov- Arnold 
intcrgals. Restricted to fl, 



7=^lnp, 



l + h 2 \q\ 2 



The level sets of / lead to all the periodic orbits (1-tori) in II. The unstable and 
stable manifolds of an 1-torus given by / = A (a constant) in A are 

W U ' S (A) = (J F u ' s {q) , 

qeA, I(q)=A 

which are three dimensional (taking into account the time dimension). 
Theorem 3.6 (Arnold Diffusion). For any A\ and A 2 such that 



— lnp <Ai<A 2 <+oo, 
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where po = 1 + h 2 (3tan-|) 2 , there exists a olq > such that when \a\ > ao, 
W u (Ai) and W S {A2) are connected by an orbit. 

PROOF. One can check directly that for any q S A, F\ (q) — — 2 and dF\ (q) /dq = 
[cf: (E3 and (|2~T2|) ]. Now consider W s {ai) and W u (a 2 ). Along any orbit q s (t) 
in W s (ai), we have 

lim 7i(g s (*)) - F!(5"(t)) = -2 - F 1 (q s (t)) 

-J-dt = -iel {F 1 ,H 1 }dt, 
lim 7(«T(t)) - 7(g = oi - /(?'(*)) 

/* +0 ° / > + °° 

= y Tt dt = -iej t {I,H l} dt, 



where 



N-i 



df_dg_ _ df dg 
dq n dq^ dcfc dq n 



is the Poisson bracket. Notice that {F\,Ho} = {I,Hq} = at any q G S. Since 
dHi/dq — > exponentially as t ^ +oo, the corresponding integrals converge. 
Similarly along any orbit q u (t) in W u (d2), we have 

F 1 (q «(*)) - lim F^q »(t)) = F x {q «(*)) + 2 

t — ► — OO 

-±dt = -ie {F 1 ,H 1 }dt, 
_ OQ at J co 

I(q u (t)) - lim 7(g »(*)) = 7(g «(*)) - a 2 

t — > — oo 

f* d7 f* 
= / -ydt = -ie \ {I,H x }dt . 

J — oo ^ J — oo 

Thus a neighborhood of W /S (a 1 ) in S can be parameterized by (7, t , t, Fx, I) where 
7 is defined in l|2.9|) . to is the initial time, and 

Fx = F 1 (q s (t))+v s 1 , 7 = 7(g s (t))+7;| . 

When yf = u| = 0, we get W s ( ai ). Thus W s ( ai ) n VF"(a 2 ) ^ if and only if 

Fxtf'®) = F 1 (q u (t)) , I(q s (t)) = I{q u {t)) , 

for some 7 and to- In such a case, there is an orbit q(t, e) C W s (ai) D W u (a2) along 
which 

/+00 r+oc 
{Fx, Hi}k(t,t)dt = , ai-a 2 = -ie {I, Hi}\^ e) dt . 
-00 J — 00 

Let <f(i, 0) be an orbit of DNLS such that q(Q, 0) and <f(0, e) have the same stable 
fiber base point. Then 

||gt0,0)-g(0,e)||~O(e) . 
For any small S > 0, there is a T > such that 



^ {F^H^t^dt <5, I* {I,Hi}\^ t) dt 

J ±00 J ±00 



< (5 , Ve S [0, e ] 
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for some eo > 0. For this T, when e is sufficiently small, 

\\q(t,0)-q(t,e)\\~O(e) , Vt € [— T, T] 

Thus 



/ {F 1 ,H 1 }\ qlt , e) dt= [ {F x ,H x }\^ m dt + 0{e) , 

J-T J-T 



Finally we have 



{Fi, ffi}|« t , e) dt = / {Fi,Hi}\xt,o)dt + 0{8) = , 

-oo J — oc 

/+oc r+OG 
{7, Hx}\^ e) dt = -ie / {I, Hx}\^ tfi) dt + 0(<5e) . 
-oc J —oc 



Next we solve the above equations at the leading order in S and e. Rewrite the 
derivative given in (|2.12|) as follows 



OFi/dqn = V n e-* , 



where 7 = 7 + 2(a 2 — ui 2 )t , t — p/fi, and V n represents the rest which does not 
depend on 7. We also rewrite q c <|2.9[) and q n (|2.11l) as 



q c = q c e 11 , q n = q„e n 



where q c and q n do not depend on 7. Then substitute all these into the leading 
order terms in H3.2fl - H3.3|) . we obtain the following equations 



(3.4) ad Ml + M| sin(t + 0i) + JM§ + M\ sin(2 7 + 2 ) = 



(3.5) ax - a 2 = 2e x / Af| + Mf sin(2 7 + 3 ) 
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where 



Mi M 2 M 3 

cos 6*i = — = , sin#i = — = , cos 6*2 = — = , 

y/M* + Ml y/Mf+Ml y/M$ + M| 

SHI #2 = — = , COS #3 = — = , SlHt/3 



VW+Ml y/Mg + Mg VW+M i 



/+00 j « j- 
^2 cost p„ Im [VnGW dr , 

,+00 N-l 

M 2 = — ^ sinr p„ Im [KG*] dr , 

J-°° n=0 

J — 00 „_n J —OO 



"°° n=0 17 -°° n=0 

M 5 = - / ^ ReG 3 n dr , M 6 = / Im G 3 n dr 

~°° n=0 ^-° c n=0 



r i _ gn+1 - 2g„ + g n _i 2 _ q n +i - 2q n + g n -i 

h 2 ' " _ /i 2 



^3 _ p gnCgn+l - 2g„ + q n -\) 
— * ^2 ' 

and r = t+p/fi. Equations (|3.4|) - (|3.5[) are easily solvable as long as neither 
y/M? + M| nor >/Mf + M| vanishes. In Figures EE we plot the graphs of them 
as functions of a. We solve equation (|3.5(l for 7, then solve equation (|3.4|) for to. 
Thus when \a\ is large enough, we have solutions. It is also clear from equations 
(I3.4|) - (|3.5|) that W s {a\) and W u (a2) intersect transversally. Then we can choose a 
sequence 

A\ = a\ < a 2 < ■ ■ ■ < ajv = A 2 , 

such that W s (aj) and W u (cij+i) (1 < j < N — 1) intersect transversally. The 
period of the 1-tori H2.H[) is ir/a 2 . Thus we can always choose the a/s such that the 
frequencies 01 the corresponding 1-tori are irrational. Therefore we obtain a 
transition chain. Apply Lemma l3.4l to the period-27r map of the DNLS, we obtain 
the claim of the theorem. □ 

Remark 3.7. The constant of motion / is equivalent to F 2 for z 2 — 1 l|2.6|l . 
The continuum limit of Hi has the form 

Hi=asint / \q x \ 2 dx + / (q 2 x + cfr 2 )dx , 
Jo Jo 

which is suitable for the NLS setting. One can regularize the perturbation by 
replacing the partial derivative d x in H± by a Fourier multiplier d x , e.g. a Galerkin 
truncation. One can use the constant of motion 



/= [ \q\ 2 dx 
Jo 



to build the second Melnikov- Arnold integral. 
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3.2. Arnold Diffusion of DNLS (N = 3, Resonant Case). In this sub- 
section, we prove the existence of Arnold diffusion for a perturbed DNLS when 
N = 3, which is a 5-dimensional system. We will study here the case that there is 
a resonance (u> > Stan^) inside the invariant plane II (|2.8|) . 
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FIGURE 3. The graph of \J Mf + M| as a function of a in the 
non-rcsonant case uj = 0. 




FIGURE 4. The graph of \/Mf + A/| as a function of a in the 
non-rcsonant case u = 0. 



Consider the following perturbation of the DNLS l|2.1|l 

H = H + e(H 1 +H 2 ) , 
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Figure 6. The graph of + M| as a function of a in the 

resonant case lo = 10. 



where 

N-l N-l 

Hi = a^(q n + <fo) , H 2 = sin t ^ 

where a is a real parameter. Under this perturbation, dynamics inside IT is changed. 
Due to the resonance a = u) in (|2.9fl . some tori do not persist into KAM tori. A 



Qn — q-n-1 
h 
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a 

Figure 7. The graph of \/Aff + Mf as a function of a in the 
resonant case w = 10. 

52 r 
50 - 
48 - 
46 - 
44 - 
42 - 
40 - 

38 I , , , , , , 

6 8 10 12 14 16 18 

a 

FIGURE 8. The graph of \J Af f + Mf as a function of a in the 
resonant case lo = 10. 

secondary separatrix is generated. Inside this separatrix are the secondary tori 
[Figure EJ. As can be seen below, resonance does not add difficulty to the Arnold 
diffusion problem. Instead of / in last subsection, we use 

H = H = H Q + eH 1 , 
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to build one of the two Melnikov-Arnold integrals. Restricted to n, the level sets 
of H produces Figure [SJ The unstable and stable manifolds of an 1-torus given by 
H = A (a constant) in A are 



W U ' S {A) 



u ^ ,s ® 

qeA, H(q)=A 



Let 



A _J_ 



2(3tan-) - — (1 + w*br) lnp 



where po = 1 + h 2 (3 tan ^) . 

Theorem 3.8 (Arnold Diffusion). For any A 1 and A 2 such that 

A* < Ai < A 2 < +oo , 

there exists a ao > such that when \a\ > ao, W U {A\) and W S (A2) are connected 
by an orbit. 

Proof. Again one can check directly that for any q £ A, F x {q) = —2 and 
dFi(q)/dq = [cf: lE3 and [(2~T2j) ]. Similar to the proof of Theorem EH consider 
W s (ai) and W u (a 2 ). Along any orbit q s (t) in W s (ai), we have 



lim F 1 (q s (t))-F 1 (q s (t)) 



Fl(q S (t)) 



dF\ 
~~dt 



dt = — i> 



/+oo 
{Fi,H 1 + H 2 }dt 



lim H(q s (t)) - H(q s (t)) = a x - H(q s (t)) 



dH 
~dt 



dt = —it 



{H,H 2 }dt 



where 



7V-1 



{f,g} =J2 Pn 



d£dg_ _ df dg 
dq n dq^ dqVi dq n 



is the Poisson bracket. Notice that {Fi,Hq} = {H,H} = at any q G S. Since 
l^t'^f ~ * ^ exponentially as t ^ +oo, the corresponding integrals converge. 
Similarly along any orbit q u (t) in W u (a 2 ), we have 

F x {q u {t)) - lim F 1 (q u (t)) = F x {q «(*)) + 2 

i— > — oo 



J —OO 



H(q u (t)) - lim iJ(«f »(<)) = H{q u {t)) - a 2 

t — > — OO 



dH 
~dt 



-dt 



{H,H 2 }dt 



Thus a neighborhood of W s (ai) in 5 can be parameterized by to, t, F Xl H) where 
$ is the angle of the 1-torus 7? = a x in IT, to is the initial time, and 



F x =F x (q s (t)) + v s x , H = H(q s (t)) + v s 2 . 
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When vf = uf = 0, we get W s (ai). Thus FF s (ai) n W u (a 2 ) ^ if and only if 

F^q^t)) = F^q^t)) , H(q s (t)) = H(q u (t)) , 

for some $ and to- In such a case, there is an orbit q(t, e) C W s (ai) n VF u (a2) along 
which 



+ OO Z' + OO 

{Fi,ffi + ff 2 }|^t,e)rft = , a 1 -a 2 = -ie {H, H 2 }\q( t ,e)dt 



Let q(t, 0) be an orbit of DNLS such that <f(0, 0) and q(0, e) have the same stable 
fiber base point. Then 

||g(0,0)-q(0,e)||~O(e) . 
For any small S > 0, there is a T > such that 



t {F u H l + H 2 }\^ e) dt <S, f {H,H 2 }\^ e) dt 



for some eo > 0. For this T, when e is sufficiently small, 



< S , Ve € [0, e ] 



||^,0) -flfce) || -0(e) , We [-T,T] 



Thus 



/ + {-Fi.-ffi + fl"a}|«(t, e )* = / + {Fx,H x + H 2 }\ sm dt + 0(e) , 

J-T J-T 

/ + {H,H 2 }\ m ^dt = [ + {H,H 2 }\ qim dt + 0(e) . 

J-T J-T 

Finally we have 

/+oo r+oo 
{F 1 ,H 1 + H 2 }\ ?(M) dt = / {F^Ht + H 2 }\ sm dt + 0(S) = , 
-oc J — oo 

/+oo /* + oo 

{H, H 2 }\ €{t ,e)dt = -it / {H , H 2 }\ ?{m dt + O(Se) 
~oo J — oo 

To the leading order terms in (|3.6)1 - H3.7|1 . we obtain the following equations 



(3.8) 
(3.9) 



\j Ml + Mf sin(t + 9i) + Q^Af| + M|sin(7 + 2 ) = 
ai - a 2 = 2e^Mf + M| sin(i + 3 ) , 
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where 



M x M 2 M 3 

— — , sin fi = — — , cos a 2 = — — 

yJM{ + Ml y/M* + Ml y/M* + M 2 



4 



sm6l 2 = —====== , cos6» 3 = —===== , sm6> 3 



y/M'i + M 2 v/Mf + M? v/Mj+Mf 



AT-l 



cost p„ Im [VnG^] dr , 
M 2 = - / ^ sinr p n Im [V n G*] dr , 

J — OC n J —CO n 



71=0 » -OO ra = Q 



N-1 



M 5 = / }J cost Im [G^G 2 ] dr , 
n =0 
/■+oo N-1 

M 6 = - sinr Im [ G » G n] dr > 

•'-OO 71=0 

^,1 _ <?Tt+i — 2g„ + g n ^i 

n — ^2 ' 

G n = T, [Qn+l - 2<ln + Qn-i I + |<7n| 2 |g n +X + <7n-i I - 2w 2 g„ , 



and r = t+p/fi. Equations (|3.8|) - (|3.9|l are easily solvable as long as neither 
yj M| + M| nor ^/ M| + M| vanishes. In Figures EE we plot the graphs of them 
as functions of a. We solve equation (|3.9|) for to, then solve equation (|3.8I) for 7. 
Thus when |a| is large enough, we have solutions. It is also clear from equations 
(I3.8|) - (|3.9|) that W s (a\) and W u (a 2 ) intersect transversally. Then we can choose a 
sequence 

Ai = a\ < a 2 < ■ ■ ■ < ajy = A 2 , 

such that W s (a,j) and W u (aj + i) (1 < j < N—l) intersect transversally. The period 
of the 1-tori (H — aj) depends on ctj no matter they are KAM tori or secondary 
tori. We can always choose the aj's such that the frequencies of the corresponding 
1-tori are irrational. We can use one secondary torus inside and close enough to 
the separatrix (Figure EJ) to bridge across the resonant region of width O(yfe). 
Therefore we obtain a transition chain. Apply Lemma l3.4l to the period-27r map of 
the DNLS, we obtain the claim of the theorem. □ 



Remark 3.9. The continuum limit of Hi and H 2 have the form 

Hi = a (q + q)dx , H 2 = suit \q x \ 2 dx , 
Jo Jo 

which is suitable for the NLS setting. One can regularize the perturbation by 
replacing the partial derivative d x in H 2 by a Fourier multiplier d x , e.g. a Galerkin 
truncation. 
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